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Abstract. We introduce a Hubbard model on a particular class of geometries, 
and consider the effect of doping the highly spin-degenerate Mott-insulating state 
with a microscopic number of holes in the extreme strong-coupling limit. The 
geometry is quite general, with pairs of atomic sites at each superlattice vertex, and 
a highly frustrated inter-atomic connectivity: the one dimensional realization is a 
chain of edge-sharing tetrahedra. The sole model parameter is the ratio of intra- 
pair to inter-pair hopping matrix elements. If the intra-pair hopping is negligible 
then introducing a microscopic number of holes results in a ferromagnetic Nagaoka 
groundstate. Conversely, if the intra-pair hopping is comparable with the inter-pair 
hopping then the groundstate is low spin with short-ranged spin correlations. We 
exactly solve the correlated motion of a pair of holes in such a state and find that, in 
1-d and 2-d, they form a bound pair on a length scale that increases with diminishing 
binding energy. This result is pertinent to the long-standing problem of hole-motion in 
the Cu02 planes of the high-temperature superconductors: we have rigorously shown 
that, on our frustrated geometry, the holes pair up and a short-ranged low-spin state 
is generated by hole motion alone. 



Submitted to: J. Phys. A: Math. Gen. 
1. Introduction 

In the years since the discovery of high-temperature superconductors, a wide spectrum 
of ideas about the fundamental cause of superconductivity in Cu02 planes has been 
proposed, but as yet none have been accepted as correct. In this article we revisit 
Anderson's early suggestion of a resonating-valence-bond (RVB) state, with 'pre-formed' 
pairs [T], which at the time lacked rigorous support. 

A suggestion currently popular in the hterature is that the three-band Hubbard 
model in the strong-coupling regime can be mapped onto the t-J model via the idea 
of Zhang- Rice singlets j2] • We believe that this reduction eliminates important physics: 
for example, the t-term alone can only generate Nagaoka ferromagnetism [3j — in order 
to account for the observed low-spin state one has to rely on the higher-order J term. 
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Figure 1. A chain of edge-sharing tetrahedra, a one-dimensional example of the class 

of geometry considered in this article. The inter-pair hopping ( ) is, in general, not 

the same strength as the intra-pair hopping ( ). 

In the accompanying article by IBS and MWL |4j, evidence is provided to support the 
hypothesis that the three-band model allows an RVB state to be generated solely by 
the hole motion. 

The three-band model for CUO2 planes is highly frustrated in a subtle way [S], 
and it is this frustration that is lost in the reduction to the t-J model. In this article 
we provide rigorous evidence for energetically favourable hole-pairing in a not dissimilar 
highly-frustrated Hubbard model in the strong-coupling limit. The frustration is crucial, 
as it is this that stabilizes the short-range correlated low-spin state, which is generated 
by hole motion alone. 

2. The model and geometry 

In this study we consider a repulsive Hubbard model on a particular class of highly- 
connected lattice. The geometry is that of a superlattice of vertex pairs. In one 
dimension we will consider a linear-chain superlattice, in two a square superlattice and 
in three a cubic superlattice. The linear chain situation, where each superlattice site 
has a co-ordination number Z = 2, can be realized physically as a chain of edge-sharing 
tetrahedra, as illustrated in figure [H 

The single-particle part of the model Hamiltonian can be specified by the options 
open to a charged spin-half fermion on any particular vertex. The fermion can hop to 
any of the vertices in the Z neighbouring pairs with a matrix element —t. It can also hop 
to the other vertex in the pair with a matrix element which, for convenience, we define 
to be —ZXt. The parameter A expresses the difference in energy scales between inter- 
and intra-pair hopping. The electrostatic correlations between the charged fermions are 
included via an on-site Hubbard term: two particles occupying the same vertex incur a 
Coulomb penalty U. The second-quantized Hamiltonian is 

^=-^J2Yl P^ra Pi'r'a ' ^ ^^Yl P^aPifa + ^ 5Z PtaPirAraPirS . (l) 

(«,«') T,T',Cr i,T,(y 

where p\^^ creates a fermion of spin a on vertex r of the pair on superlattice site i. The 
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complementary spin, or vertex, is denoted by a bar, and {i, i') indicates summation over 
nearest neighbours i and i'. 

We are concerned with behaviour adjacent to the Mott-insulator, and so it is 
appropriate to rewrite the Hamihonian in terms of fermionic hole operators, via the 
mapping p]^^ ^jrcr- Then, up to an ignorable constant, the hole Hamiltonian is 

^ = ^ E E h.r'a + ZXt J2 hlaKa + ^ E hLKrAr.Kr. ■ (2) 

In this article we restrict attention to the extreme strong-coupling limit of f/ = cxd 
with a finite positive t. At maximal filling in the finite energy subspace the system is 
a Mott insulator with one particle (or, cquivalcntly, one hole) per vertex: this state is 
hugely spin-degenerate because the super-exchange energy scale, J, is rigorously equal 
to zero. We address the problem of how this degeneracy is lifted by the inclusion of a 
microscopic number of holes. This particular geometry is useful because we can treat 
the doping of one and two holes analytically. Before the more interesting issue of a 
two-hole bound state we will introduce the two competing styles of groundstate via the 
more straightforward single-hole situation. 

3. Local symmetries and one-hole behaviour 

Eigenstates respect the symmetries of the Hamiltonian, and this is useful when 
categorizing the possible styles of solution. For example, the model that we consider 
in this article preserves the spin of a hole when it hops: as a result of this global 
symmetry, the total spin of a state is conserved. There also exists a collection of local 
symmetries, where N is the number of superlattice sites: these arise from the fact that 
the Hamiltonian is unchanged when the two vertices of any pair are swapped. In looking 
for the groundstate of our Hamiltonian we categorize the eigenstates by their symmetry 
with respect to these N Hamiltonian-preserving operations. We consider states that are 
either symmetric or antisymmetric with respect to all N of the local symmetries: we do 
not deal with the mixed-symmetry states. 

3.1. The high-spin state 

Setting the parameter A to zero effectively removes all the intra-pair bonds; if the 
superlattice is bi-partite then the connectivity is no longer frustrated, and hence the 
one-hole groundstate will exhibit Nagaoka ferromagnetism. 

In the real-space representation that we choose to work with, every site is occupied 
with either a spinless fermionic hole or a hard-core bosonic spin. To illustrate this, 
consider a small part of the geometry comprising two nearest-neighbour vertex pairs in 
the ferromagnetic one-hole state: if the hole is on one of the vertices then the other three 
vertices accommodate parallel spins. In pictorial form, the action of the Hamiltonian 
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on this unit is simply 

h — ^er 



H 



ZXt 



a a 



a 



+ t 



a, — a- 
+ 



(3) 



We can use this representation to write down the Bloch state for one hole with 
momentum A; in a ferromagnetic background: 




(4) 

-1 i i+1 i— 1 i i+1 

The one-dimensional lattice is used for illustration here because it is simple to depict; 
it is obvious how this representation extends to the two- and three-dimensional lattices. 
This state is 'fully symmetric' since it maps onto itself when the vertices of any particular 
pair are interchanged. 

With the aid of equation it is straightforward to show that \Pko) is an eigenstate 
of the Hamiltonian with an energy eigenvalue 

,(higM)^^ + 27.. (5) 

Throughout this article, energy eigenvalues, denoted by e, are 'per hole' and in units 
of Zt. The one-hole groundstate has the wavevector k that minimizes the superlattice 
structure factor 7^: for the cubic family of superlattices that we are considering here 
this means that fc = Q = (tt, tt, . . .) and therefore 



^(high,i) 

Qcr 



A-2. 



(6) 



3.2. The low spin state 

The other class of state that we consider is the complement of the high-spin state, in that 
it has low-spin and is 'fully anti-symmetric' with respect to the vertex-pair-swapping 
symmetry. 

In our chosen representation any real-space pair that is free from holes must 
accommodate two bosonic spins: in order to make this pair of bosons antisymmetric, the 
spins must be correlated in a total-spin singlet. We depict the singlet in the following 
way: 



v/2 



(7) 



The antisymmetric property of the bosonic singlet is evident from the fact that swapping 
the 'head' and 'tail' generates a factor of —1. 

The one-hole antisymmetric low-spin state is 



|Lfco-) 



.ifc.ri 



1 



V2 





(8) 



i— 1 i i+1 i— 1 i i+1 

where the anti-phase combination guarantees antisymmetry with respect to vertex 
exchange on the hole's pair. 
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The action of the Hamiltonian on a two-pair unit that contains the hole is 



H 



h — 71 a I r — ^h 



a h 



(9) 



Acting the Hamiltonian on the state |Lfco-) creates states where singlets extend from one 
superlattice site to another. Such states have a non-zero overlap with the original state, 
and we can use the spin identity 



to show that |Lfccr) is in fact an eigenstate of the Hamiltonian with energy eigenvalue 

-A + 7fc. (11) 



(low,l) 



As was the case for the ferromagnetic one-hole state, the lowest energy state has k = Q. 
By comparing equations © and (fTTj) it is obvious that the groundstate is ferromagnetic 
when A < 1/2 and low-spin with short-range correlations when A > 1/2. 



4. Two-hole behaviour 



In the previous section we saw that the one-hole solutions were trivial to find once a 
suitable representation had been established. With this framework in place, we are now 
ready to tackle the two-hole problem. The high-spin state is simple to describe because 
the interaction is irrelevant and we simply occupy a single non-interacting band; the 
low-spin problem is not so straightforward, and we make use of the Greens function 
impurity technique to deal with the interaction. 



4-.1. The high spin state 

The single hole problem in the symmetric subspace was solved by including a symmetric 
hole with wavevector Q; all other sites feature a spin a hence the state was described as 
'high spin'. If we include another symmetric hole from the corner of the Brillouin zone, 
with the same spin as the first, then Pauli exclusion prevents double-occupancy of any 
vertex. On a macroscopic lattice these two holes, which are in orthogonal single-hole 
states, can have an arbitrarily close energy. Hence the lowest energy per hole of the 
two-hole high-spin state is the same as the single-hole minimum, namely 



4-2. The low spin state 

Here we restrict attention to states where the pair of fermionic holes, hi and h.2, have 
zero total momentum. The other degree of freedom is the separation between the two 
holes on the superlattice, R = r2 — vi, and we use this quantity to label the low-spin 
state basis. 
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Consider the R = state which, in its properly normahsed form, is 



\R = 0) 




(13) 



^ i-1 i i+1 

This state has total-spin zero, short-range spin correlations, and is overall antisymmetric 
since swapping a singlet's head and tail, or exchanging hi and h.2, result in a minus sign. 

As before, consider a unit of this state comprising the vertex pair containing hi and 
h.2, and one of its nearest neighbour pairs that must feature a singlet. The action of the 
Hamiltonian on this unit can be written as 



h^ 
li^A 



t 



-ill 



li^> 



+ 



h2-^i 



+ 



(14) 



oo 



Note that no use is made of the intra-pair hopping, as this is prohibited by the U 
Coulomb penalty. 

Now, any particular superlattice site has Z nearest neighbours; let us denote the 
Z vectors that translate to the nearest neighbours as ni, n2, . . . , nz- We then choose 
to define the overall antisymmetric two-hole state where the holes are separated by one 
lattice spacing as 

" T h2 hi i 1— i _ 



h2 hi i 
hi i ^ t ii2 




(15) 



where it is implied that all other vertex-pairs feature a spin singlet. The convention is 
to put an I spin on the same pair as the hole hi, and a | with h2; the factor of 1/2 
ensures correct normalisation. It is then possible to express the outcome of applying 
the Hamiltonian to |-R = 0) as 

z 

H\R = 0) = 2^y2t'^\R = nj). (16) 



The states for all other R, i.e. R ^ 0,ni, n2, . . . , nz, are defined in a similar way 
to \R = rij) . The fact that these states have the holes more than one hop apart means 
that they move as 'independent holes'. The reasoning from section applies, and thus 



z 



H\R^0,ni,n2,...,nz) = -2ZXt\R) + 2ty \R + nj) 



(17) 



where the factor of two arises because there are two independent holes. 

The action of the Hamiltonian on the nearest-neighbour states is a combination of 
the two above styles of behaviour. The holes are on different vertex pairs, so they can 
gain the intra-pair hopping energy. Hopping in one particular direction will bring both 
holes onto the same site, and hopping in all the other directions is 'independent-hole'- 
like. The result is that 



H\R = np) = -2ZXt\R) + 2y/2t\0) + 2t^\R + nj) . 



(18) 
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Now that we know all the matrix elements, it can be seen that they deviate from 
a perfectly periodic Hamiltonian in just two ways: (i) the coupling between |0) and 
the 'nearest-neighbour' states \nj) is enhanced by a factor of y/2, and (ii) the intra-pair 
energy 2ZXt is not gained by |0). We can write the Hamiltonian as the sum of a periodic 
term, Hq, and a 'local' term Hi, where the matrix elements of these two terms are given 
by 

z 

Ho\R) = -2Z\t\R) + 2t'^\R- 



(19) 



and 





' |o) ■ 




' z\ 72-1 72-1 ■ 


■ 72 - 1 " 




" |o) ■ 




l^i) 




72-1 







1^1) 


Hi 




= 2t 


72-1 











. \i^z) . 




.72-1 










(20) 



It is convenient to rescale the Hamiltonian so that all energies are 'per hole' in 
units of Zt: the rescaled two-hole Hamiltonian, Ti, is defined hj H = 2ZtT-C. In order 
to exactly solve this eigenproblem we will make use of the Greens function for the entire 
Hamiltonian, 7i = Tig + Ti-i, which can be written as 

g{z) = go{z) + go{z)j:{z)go{z) (21) 

where go{z) is the Greens function for T-Cq, the periodic part of the Hamiltonian, 

g,{z) = [z-no]-' = ^ (22) 

^ + A - 7fc 

and 

E{z) = [i-nigo{z)]-'ni. 



(23) 

The energy eigenvalues of the periodic case alone are found at the poles of go{z). If we 
anticipate that the inclusion of Hi generates a new state, then its pole must be present 
in E(z). Hence the pole is found at z = e such that 

det [1 - Higoie)] = 0. (24) 

Equation (|2(J|) shows how the local term of the Hamiltonian, Tii, can be represented by 
a (Z + 1) X (Z + 1) matrix; the above determinant equation involves the product Tii^O) 
hence the relevant part of the matrix go{z) is simply 

^00(2;) goi{z) ^02 (^) 
^10(2;) ^11(2;) guiz) 

go{z)= Q20{Z) g2l{z) g22{z) ■■■ g2z{z) (25) 

. gzo{z) gzi{z) gz2{z) 



■ ■ ■ g^z 


z) 1 


■■■ giz 




■ ■ ■ g2z 




■ ■ ■ gzz 
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where, if we conveniently let no = 0, these real-space elements are defined by 

1 exp[i A;.(nj — rij)] 



Gijiz) 



Extracting the point-symmetry of the superlattice reduces the matrices to 2x2, 



where 



Hi 



A 

v/2-1 



v/2- 1 



N ^ 



Ik 

2; + A - 7fc 



and Qq 



and Q{z) 



-E 



2; + A - 7fc 



Evaluating the determinant equation (|24j) to find the pole at 2; = e gives 

1 = A^oo(e) + 2 (72 - 1) Q\t) + (^2 - 1)^ [^oo(e) Q[e) - Q'\e)\ 

which can be further simplified using the identities 

6?'(e) = (e + A)6?oo(e)-l and Q{e)G^^{e) 

to give the self consistent equation 

N ^ 



1 



e + 2A V e + A - 7fc 



(26) 

(27) 
(28) 

(29) 
(30) 

(31) 



In one- and two- dimensions this equation has a solution in the parameter range 
< A < 1. For example, the energy per hole of the linear-chain superlattice state is 
given exactly by 

2 r, rr-. rl _^ _ ^ _ ^ ^ ^32) 



a + Ta^Ts 



where we have parameterized the binding energy per hole as 5: this energy is compared 
with the independent-hole low-spin and high-spin energies in figure El Equation 1)311) 
can be solved numerically for the case of the square superlattice; the binding energy is 
about an order of magnitude smaller than the 1-d binding energy, and hence the 2-d 
bound state curve is barely visible below the independent-hole fine in figure El 

The Greens function contains the eigenvalue and eigenvector information, so this 
technique allows us to determine the 'size' of the bound-state pair as a function of 
binding energy. Provided that 5 7^ 0, the amplitude of the real-space wavefunction 
when the holes are separated by Ri is proportional to 

l_ ^ exp[ifc.J^i] ^^^^ 



^io(e) 



e + A - 7fc 



In the one-dimensional example, this quantity is 



3 (A + yA"2 



(34) 



2A- v^FTS 

which means that, for A < 1, the magnitude of the real-space wavefunction decays 
exponentially with a length scale ^ defined by 
1 



In 



i(A + v^ 



(35) 



The variation of this length scale with A is plotted in figure (21 
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Figure 2. The energy per hole of various states as a function of the model parameter 
A. Whether we are considering the linear chain, square, or cubic superlattice, the 

one-hole ferromagnetic Nagaoka state ( — • — ) and one-hole low-spin state ( ) are 

degenerate at Ac^^ = 1/2. The two-hole bound-pair low-spin state on the linear chain 

( ) is degenerate with the two-hole Nagaoka state at Ac — 0.4809(6). The binding 

energy of the two-hole low-spin state on the square superlattice ( ) is an order of 

magnitude smaller than the linear chain case; in both 1-d and 2-d the binding energy 
vanishes at A = 1. 

5. Discussion and conclusions 

We have shown that, for a restricted range of A, pairs of holes in our strong-coupling 
Hubbard model bind together to form the real-space analogue of Cooper pairs. The 
background spins are induced into an RVB state by the motion of these holes. We have 
not addressed the important issue of the relevance of these results to superconductivity 
in Cu02 planes, and whether or not one would anticipate superconductivity in our model 
at finite doping 

It is important to realize that the demands of superconductivity are less stringent 
than that of binding holes into pairs. In one and two dimensions any magnitude of net 
attraction between two holes guarantees a bound-state, but in three dimensions there is 
a non-zero threshold that must be exceeded. Consequently, we can achieve the bound- 
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Figure 3. The size of the bound-state hole pair ^ (in units of the lattice parameter) 
on the one- dimensional chain, as a function of the model parameter, A. The binding 
energy vanishes as A — > 1, with a concomitant divergence in the real-space size of the 
bound-pair. 

state for the linear chain and square superlattice, but in the cubic superlattice case the 
bound-state is irrelevant as it only forms for values of A for which the groundstate is 
ferromagnetic. 

Superconductivity, on the other hand, involves macroscopic doping. When the 
number of holes is of the same order as the number of lattice sites, the holes are forced 
to be near each other, and then can always make use of an attraction to form a 'collective' 
paired state: attraction guarantees a pairing instability of a Fermi surface. It is possible 
that in one dimension a Luttinger liquid might arise, with power-law pairing correlations 
and no long-range order, but inter-layer coupling stabilizes the superconductor. 
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